Tunable Superconducting Phase Transition in Metal-Decorated Graphene Sheets 
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Using typical experimental techniques it is difficult to separate the effects of carrier density from 
disorder on a two-dimensional superconducting transition. To address this problem, we have pro- 
duced graphene sheets decorated with a non-percolating network of nanoscale tin clusters. These 
metal clusters both efficiently dope the graphene substrate and induce long-range superconducting 
correlations. This allows us to study the superconducting transition at fixed disorder and vari- 
able carrier concentration via the field effect. Wc find that despite structural inhomogeneity on 
mesoscopic length scales (10-100 nm), this material behaves electronically as a homogenous dirty 
superconductor. Our simple self-assembly method establishes graphene as an ideal tunable sub- 
strate for studying induced two-dimensional electronic systems at fixed disorder and our technique 
can readily be extended to other order parameters such as magnetism. 



The superconducting transition in two-dimensions is 
of interest for both the fundamental understanding of 
electronic order in reduced dimensions and applications 
involving superconducting thin films. An open question 
is how the transition behaves as a the density of carri- 
ers mediating the superconductivity is varied. In par- 
ticular, the strength of disorder appears to play a funda- 
mental role separating qualitatively different behavior[T]. 
In two-dimensions, the electric field-effect provides the 
most versatile method for tuning the carrier density of 
a system at fixed disorder. However, the field effect 
places more stringent limits on the dimensionality of the 
system since the film must be thinner than the Debye 
length, which is typically much smaller than the pene- 
tration depth that places the limit on two-dimensional 
superconductors. Despite this limitation, the field effect 
has been used to tune the superconducting transition in 
specific materials such as thin films with anomalously low 
carrier density [2] and interfacial states between complex 
oxides [3]. An alternative approach is to couple supercon- 
ducting correlations directly into a truly two-dimensional 
electronic system, graphene. 

The bipolar two-dimensional electron gas (2DEG) 
present in graphene |4j is markedly different from the 
buried 2DEGs found at oxide interfaces or in GaAs het- 
erostructures in that it is 'open' to the environment with 
a stable and inert surface. Using standard 2DEGs it 
is only feasible to capacitively couple the electron gas 
to materials deposited on their surface [5) |6], whereas 
graphene is expected to allow direct coupling and thus 
offer access to different regions of phase space [7]. We 
sought to determine whether bulk materials deposited 
directly onto the graphene surface, such as metal clus- 
ters, can act as dopants and efficiently couple through 
the electron gas, whose carrier density and type can be 
tuned by an applied gate voltage. The low carrier den- 
sity in graphene, relative to bulk values, and weak in- 



trinsic interactions such as spin-orbit coupling, should 
limit backaction of the electron gas on dopant proper- 
ties. Thus, exposed graphene sheets could provide a near- 
ideal substrate for the manipulation and general study of 
proximity-induced electronic phases. 

Graphene has been shown to effectively carry 
proximity-induced Josephson currents injected from con- 
tacting electrodes [HI However, a finite coherence 
length limits the length of such junctions to approxi- 
mately one micron, reducing the physics to one dimen- 
sion. To maintain coherence over longer distances in two 
dimensions while retaining the unique properties of the 
graphene sheet we employ a geometry (Fig. la) where 
a large array of nanoscale dopant islands is placed in a 
non-percolating network on top of the graphene sheet [TU]. 

We avoid complicated lithographic patterning and ex- 
ploit the poor wettability of graphite to simply and reli- 
ably produce an array of submicron islands. Low melt- 
ing point metals such as the elemental superconductor 
Sn readily form self-assembled islands when deposited on 
pristine graphene at room temperature (Fig. lb) similar 
to previous results on graphite ^JJ. Analysis of scanning 
electron micrographs and atomic force micrographs in- 
dicates that 10 nm of nominal deposition thickness typi- 
cally results in islands with 80±5 nm diameter and 25±10 
nm gaps between them (Fig lb). In general, many dif- 
ferent materials with different electronic order parame- 
ters can be deposited via this process by controlling the 
graphene substrate temperature during deposition |12j. 
and other deposition methods such as chemical function- 
alization and wet self-assembly could be used as well. 

Samples were prepared by exfoliating Kish graphite H] 
onto degenerately doped {p < 0.005 mfi-cm) silicon 
wafers coated with 285 nm of thermal oxide. Single- 
layer graphene flakes were identified by optical con- 
trast and confirmed via analysis with a micro-Raman 
spectrometer [TB] . Four-probe contacts were defined via 
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FIG. 1: (Color online) a) Schematic of device configuration 
and measurement setup. Blue islands correspond to Sn clus- 
ters, b) Scanning electron micrograph of Sn island morphol- 
ogy on the graphene sheet (Scale bar — 100 nm. Inset: optical 
image of a typical device showing the four probe configura- 
tion (Scale bar 10 microns) c) Four-terminal sheet resistance 
as a function of gate voltage for Sample A before (red online) 
and after (blue online) Sn deposition. The dotted lines in- 
dicate the charge neutral point and the arrow indicates the 
shift after Sn deposition. 



electron beam lithography and a Pd/Au 10/50 nm bilayer 
was evaporated as metal electrodes (Fig. 1, c). To pro- 
duce the island network, Sn (99.999% purity) was evap- 
orated using an electron gun in high vacuum (10~^ torr) 
onto graphene substrates at room temperature. The sam- 
ples were thermally anchored to the cold stage of He-3 or 
Dilution cryostats and connected to highly filtered lines. 
Linear response and differential conductance were mea- 
sured with standard low frequency lock-in techniques us- 
ing low excitation currents in the range 10-100 nA. Figure 
Ic displays the room-temperature field-effect characteris- 
tics of a device before and after Sn is deposited. Although 
40% of the graphene surface is coated by Sn islands after 
the deposition, many of the original electronic properties 
of graphene remain intact, including bipolar transport 
and field-effect mobilities fi > 1000 cm2/(V • s). The 
three main effects of Sn deposition are a rigid shift in 
the charge neutrality point (Dirac point Vd) to more 
negative voltages, a factor of five decrease in mobility, 
and a pronounced asymmetry between electron and hole 
transport. Note that the maximum resistance at charge 
neutrality remains unchanged, indicating that the high 
coverage of low resistance Sn islands does not directly 
shunt the current in the graphene sheet. 

All three of these effects are well described by inho- 
mogeneous doping due to charge transfer from the metal 
islands to the graphene sheet [14]. From the shift in the 
charge neutrality point (Vd) on three separate samples 
and the known gate capacitance {Cg = 115 aF / /xm^) we 
can calculate the charge induced in the graphene sheet 



by the Sn, Uind = CgAVo- Normalizing by the observed 
Sn coverage, we infer that Sn transfers 9 ± 2 x 10^^ 
cm^^ electrons to the graphene underneath it. This is 
expected from the difference in work functions between 
the two materials ($g = 4.5 eV, <i>s„ = 4.42 eV)[TS] 
and in agreement with recent experiments performed us- 
ing other metals [H]. This induced charge reduces the 
mobility of both types of carriers via charged impurity 
scattering [T7] while the asymmetry in transport occurs 
because holes experience the pinned Fermi level under 
the Sn islands as a potential barrier, while electrons ex- 
perience a potential well[TU[TS]. 

More interesting than the infiuence of the Sn islands 
on the normal state properties of graphene is the effect 
the superconducting correlations in the Sn have on trans- 
port via the proximity effect. Figure 2 shows the sheet 
resistance versus temperature for gate voltages on both 
the hole (Fig. 2a) and electron (Fig. 2b) sides of the 
charge neutrality point. Each curve exhibits two distinct 
features, a high temperature partial drop in resistance 
that occurs at ^ 3.5 K independent of gate voltage, and 
a broad transition between 3 K and 1 K to a state of zero 
resistance that is strongly dependent on the gate voltage. 

The first partial resistance drop (arrows in Fig. 2a, b) 
is due to condensation of Cooper pairs in the Sn islands 
(^rj^Buik ^ 3 K). Analysis of the drop shows that it 
can be fit by Aa{T) cx ln(T/rco)"^ (Fig. 3a) typical of 
Aslamazov-Larkin fiuctuation-enhanced conductivity [T5] 
in two dimensions. Fitting each curve from 3.8 — 4.5 
K we extract a mean-field pairing temperature (Tco) of 
3.54 ± 0.02 K independent of gate voltage (solid squares 
in Fig. 3c). Note that the amphtude of this drop is not 
directly proportional to the Sn coverage, indicating that 
the islands do not act as simple superconducting shunts. 
In two-dimensional superconducting systems, it is well 
known that although the amplitude of the superconduct- 
ing wave function is well defined below the pairing tem- 
perature Tco, thermally induced phase fiuctuations (vor- 
tices) destroy global phase coherence and produce dissi- 
pation due to a finite fiux fiow resistance [201 HI]' How- 
ever, below the critical Berezinskii-Kosterlitz-Thouless 
unbinding temperature, Tbkt, the attractive interaction 
between vortices with opposite orientation causes them 
to form bound pairs allowing a finite supercurrent to fiow. 

The vortex-unbinding temperature can be identified 
from the universal form of the flux flow resistance [22] 
above the transition Ru{T) cx exp[6(T — Tbkt)~^^^], 
where 6 is a constant of order unity governing the vortex- 
antivortex interaction strength and Tbkt is the vortex 
unbinding temperature. To extract this form we plot 
(dln(-Rn)/(iT) "^^^ vs. T, which produces a straight line 
with Tbkt given by the x-intercept for curves follow- 
ing the universal form. In Fig. 3b the resulting fits are 
given showing Tbkt extracted from the x-intercepts for 
three different gate voltages. Through a large intermedi- 
ate range the curves follow the universal form (straight 
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FIG. 2; Sheet resistance versus temperature for various 
gate voltages, Vg, referenced to the charge neutrahty point 
Vd = +40 V for this device. In a) AVg ^ Vg - Vd < 
corresponds to hole transport, whereas, AVg > in b) cor- 
responds to electron transport through the graphene sheet. 
The arrow labeled Tco indicates the first partial resistance 
drop corresponding to the mean-field pairing transition of the 
Sn islands. Inset: Same data on a log scale. 



line). However, at low temperatures the curves level 
off due to finite-size effects, which cut off the attractive 
vortex-antivortex interaction [22 . This departure from 
the universal form is particularly evident near the charge 
neutrality point (Fig. 2a and b), possibly indicating prox- 
imity to a superconductor-normal quantum critical point 
at the lowest charge densities [23] . Figure 3c summarizes 
the resulting Tbkt (open circles) extracted from the re- 
sistance versus temperature curves at each gate voltage. 

To analyze the gate voltage dependence of the vortex 
unbinding transition, we consider the sheet as a dirty 
two-dimensional superconductor where the gate voltage 
allows us to tune the normal state resistance. This is 
justified since the length scales of disorder are much 
smaller than the superconducting coherence length, ie. 
r ~ c? ~ Imfp < Co I where r is the size of the islands, 
d is the distance between islands, imfp '--^ 20 — 30 nm is 
the mean free path extracted from field effect measure- 
ments at 6 K and the superconducting coherence length 
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FIG. 3: (Color online) a) Fits of the sheet resistance ver- 
sus temperature to fluctuation-enhanced conductivity of the 
Aslamazov-Larkin form, b) Rescaling of the sheet resistance 
versus temperature to the BKT form to extract the vortex 
unbinding temperature Tbkt c) The mean-field pairing tem- 
perature, Tco (black squares) and vortex-unbinding tempera- 
ture, Tbkt (open circles) as a function of gate voltage. The 
solid line is a fit of Tbkt using equation f and the measured 
normal state properties of the device, (see text). 



Sn (Cq" ^ 300 nm). For a dirty 2D superconductor, one 
can use the jump in superfluid stiffness at the vortex un- 
binding transition to relate Tbkt to the normal state 
resistance of the film |24l 125] . 
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where A(r) is the superconducting energy gap, Rn is 
the normal state sheet resistance, Rq — ^•"'"f ^ « 8.96kQ, 
and €v is an effective dielectric constant that describes 
the material dependent screening of the attractive vortex- 
antivortex interaction [25 . Using the weak-coupling BCS 
limit for the superconducting gap and the sheet resistance 
measured at 6 K, we fit the Tbkt extracted above using 
ey as the only adjustable parameter and find = 2 AO ± 
0.05 (solid curve in Fig. 3), which is twice the value seen 
in related systems [25]. This implies that the vortices and 
antivortices are relatively weakly bound in our system. 

At temperatures below Tbkt, vortices and antivor- 
tices form bound pairs and a finite critical current de- 
velops which saturates to a gate- voltage dependent value 
for T <C Tbkt- Current- voltage characteristics at 100 
niK for different applied gate voltages are shown in Fig- 
ure 4. The gate-tunable critical current is qualitatively 
similar to isolated graphene Josephson junctions[Hl [S] 
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FIG. 4: (Color online) Current-voltage curves taken at zero 
magnetic field and 100 mK, corresponding to hole transport 
in the graphene sheet for gate voltages, AVg relative to the 
Dirac point Vd = +12 V for this sample. Inset: similiar 
curves for electron transport. 



with the exception that, in our devices, critical current 
densities {Ic/width 1 A/m) comparable to submicron 
graphene Josephson junctions are maintained over dis- 
tances of tens of microns, demonstrating the fully two- 
dimensional phase coherence in this system. 

We have demonstrated a simple method to produce a 
two-dimensional superconductor on a graphene substrate 
and tune the transition via an electrostatic gate. This 
allowed us to systematically tune the carrier density at 
fixed disorder. Although structurally inhomogenous, this 
material behaved electronically as a weakly disordered 
two-dimensional superconductor. While we have probed 
the properties of this system using electron transport, the 
readily accessible interface allows application of a myr- 
iad of local characterization techniques such as scanning 
probe microscopy, optical spectroscopy, etc. We expect 
arranging the islands into regular arrays or superlattices 
should lead to interesting frustration effects as a func- 
tion of applied magnetic and electric fields |26l [27] . This 
tunable superconducting material may find applications 
in bolometers for sensing applications or as an element 
in circuits for quantum information processing. The ease 
of fabrication and considerable versatility of deposition 
materials make graphene an attractive platform for in- 
vestigating other electronic orders such as magnetism in 
two dimensions. 
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